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Abstract — We consider the 2x2x2 multiple-input multiple- 
output interference channel where two source-destination pairs 
wish to communicate with the aid of two intermediate relays. 
In this paper, we propose a novel lattice strategy called Aligned 
Precoded Compute-and-Forward (PCoF). This scheme consists 
of two phases: 1) Using the CoF framework based on signal 
alignment we transform the Gaussian network into a deterministic 
finite field network. 2) Using linear precoding (over finite field) we 
eliminate the end-to-end interference in the finite field domain. 
Further, we exploit the algebraic structure of lattices to enhance 
the performance at finite SNR, such that beyond a degree of 
freedom result (also achievable by other means). We can also 
show that Aligned PCoF outperforms time-sharing in a range of 
reasonably moderate SNR, with increasing gain as SNR increases. 



I. Introduction 

In recent years, significant progress has been made on the 
understanding of the theoretical limits of wireless commu- 
nication networks. In JT), the capacity of multiple multicast 
network (where every destination desires all messages) is 
approximated within a constant gap independent of SNR and 
of the realization of the channel coefficients. Also, for multiple 
flows over a single hop, new capacity approximations were 
obtained in the form of degrees of freedom (DoF), generalized 
degrees of freedom (GDoF), and O(l) approximations Q-Q. 
Yet, the study of multiple flows over multiple hops remains 
largely unsolved. The 2x2x2 Gaussian interference channel 
(IC) has received much attention recently, being one of the 
fundamental building blocks to characterize the DoFs of two- 
flows networks [5| One natural approach is to consider this 
model as a cascade of two ICs. In JSJ, the authors apply 
the Han-Kobayashi scheme [7] for the first hop to split each 
message into private and common parts. Relays can cooperate 
using the shared information (i.e., common messages) for 
the second hop, in order to improve the data rates. This 
approach is known to be highly suboptimal at high signal- 
to-noise ratios (SNRs), since two-user IC can only achieve 
1 DoF. In ||8j, Cadembe and Jafar show that | DoF is 
achievable by viewing each hop as an X-channel. This is 
accomplished using the interference alignment scheme for 
each hop. Recently, the optimal DoF was obtained in [9] using 
aligned interference neutralization, which appropriately com- 
bines interference alignment and interference neutralization. 




Fig. 1. 2x2x2 MIMO Gaussian interference channel. 



Also, there was the recent extension to the K x K x K 



Gaussian IC in |10|, achieving the optimal K DoF using 
aligned network diagonalization. 

We consider the 2x2x2 multiple-input multiple-output 
(MIMO) IC as shown in Fig. [T] consisting of two sources, 
two relays, and two destinations. All nodes have M multiple 
antennas. Each source k has a message for its intended 
destination fc, for k = 1,2. In the first hop, a block of 
n channel uses of the discrete-time complex MIMO IC is 
described by 



(1) 



where the matrices X fe Q and Y Rfc contain, arranged by rows, 
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the fc-th source channel input sequences x fc g G 



for £ = 



1 , . . . , M, the fc-th relay channel output sequences y 



ilxn 



and where F,fc G 



iMxM 



denotes the channel matrix 



from source k to j. In the second hop, a block of n channel 
uses of the discrete-time complex MIMO IC is described by 
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It and Y fe contain the fc-th relay channel 
xn , the fc-th destination channel 

lxn , , jiMxM 



where the matrices X 
input sequences x Rfc { G 

output sequences y fc G C lx ", and where Gkj G 
denotes the channel matrix between relay j and destination fc. 
The matrix Z k (or Z Rfc ) contains i.i.d. Gaussian noise samples 

1 We use "underline" to denote the matrices whose horizontal dimension 
(column index) denotes "time" and vertical dimension (row index) runs across 
the antennas. For any matrix X, we denote x f by the £-th row of X. 



~ CM (0,1). We assume that the elements of Fj^ and Gkj 
are drawn i.i.d. according to a continuous distribution (i.e., 
Gaussian distribution). The channel matrices are assumed to 
be constant over the whole block of length n and known to 
all nodes. Also, we consider a sum-power constraint equal to 
A/SNR at each transmitter. 

In this paper, we propose a novel lattice strategy named 
Aligned Precoded Compute-and-Forward (PCoF). CoF makes 
use of lattice codes so that each receiver can reliably decode a 
linear combination with integer coefficients of the interfering 
codewords. Thanks to the fact that lattices are modules over 
the ring of integers, the linear combinations translates directly 
into a linear combination of messages over a suitable finite 
field. In this way, each hop is transformed into a deterministic 
noiseless finite field IC. The end-to-end interferences in the 
finite field domain are eliminated by distributed precoding 
(over finite field) at relays. Using this framework, we char- 
acterize a symmetric sum rate and prove that 2M — 1 DoF 
is achievable by lattice coding (this DoF result is proven in 
[9 1 without resorting to CoF and lattice coding). Further, we 
use the lattice codes algebraic structure in order to obtain 
also good performance at finite SNRs. We use integer-forcing 
receiver (IFR) of fTT) in order to minimize the impact of noise 
boosting at the receivers, and integer-forcing beamforming 
(IFB), proposed by the authors in [12) , in order to minimize the 
power penalty at the transmitters. We provide numerical results 
showing that Aligned PCoF outperforms time-sharing even at 
reasonably moderate SNR, with increasing performance gain 
as SNR increases. 

II. Preliminaries 

In this section we provide some basic definitions and results 
that will be extensively used in the sequel. 

A. Nested Lattice Codes 

Let Z[j] be the ring of Gaussian integers and p be a prime. 
Let © denote the addition over ¥ q with q = p 2 , and let g : 
F q — > C be the natural mapping of ¥ q onto {a + jb : a,b G 
Z p } C C. We recall the nested lattice code construction given 
in (T3). Let A = {A = zT : z G Z n \j]} be a lattice in C\ 
with full-rank generator matrix T G C™ x ". Let C = {c = 
wG : w e FJ} denote a linear code over ¥ q with block length 
n and dimension r, with generator matrix G. The lattice Ai 
is defined through "construction A" (see fl4) and references 
therein) as 

A 1 =p- 1 g(C)T + A, (3) 

where g(C) is the image of C under the mapping g (applied 
component-wise). It follows that A C Ai C p~ 1 h is a chain 
of nested lattices, such that |Ai/A| = p 2r and |p _1 A/Ai| = 

p 2(n-r)_ 

For a lattice A and r 6 C , we define the lattice quantizer 
Qa(e) = argmin^ eA ||r — A|| 2 , the Voronoi region Va = 
{r G C n : Q A (r) = 0} and [r] mod A = r - Q A (r). For A 
and Ai given above, we define the lattice code C = Ai n Va 
with rate R = -log|£| = ^logg. Construction A provides 



a natural labeling of the codewords of L by the information 
messages w e ¥ r q . Notice that the set p~ x g(C)T is a system 
of coset representatives of the cosets of A in Ai. Hence, the 
natural labeling function / : — >• C is defined by /(w) = 
P~\g(wG)T mod A. 

B. Compute-and-Forward 

We recall here the CoF scheme of fPJ| . Consider a 2-user 
Gaussian multiple access channel (MAC) with M antennas at 
each transmitter and at the receiver, represented by 

Y = HCX + Z (4) 

where H G C MxM , C G Z[j] Mx2M , and Z contains i.i.d. 
Gaussian noise samples ~ CN(0, 1). This particular form of 
channel matrix HC will be widely considered in this paper, 
as a consequence of signal alignment. All users make use 
of the same nested lattice codebook C = Ai n Va, where 
A has second moment o\ — y ^ y : J Vx ||r|| 2 <ir = SNR cfr . 
Each user k encodes M information messages w fe t G ¥ r q into 
the corresponding codeword t ke = /(w fc e ) and produces its 
channel input according to 

*k,e = [h,e + &kA m °d A, (5) 

for i = 1,...,M, where the dithering sequences d fe £ 's 
mutually independent across the users and the messages, 
uniformly distributed over Va, and known to the receiver. 

Notice that Xfc i i s me row °f 2^fc' ^ = 1)2. The 
decoder's goal is to recover an integer linear combination 
s = [b H CT] mod A, with some integer vector b € Z[?] Mxl . 
Since Ai is a Z[j]-module (closed under liner combinations 
with Gaussian integer coefficients), then s G C. Letting s 
be decoded codeword, we say that a computation rate R is 
achievable for this setting if there exists sequences of lattice 
codes C of rate R and increasing block length n, such that the 
decoding error probability satisfies lira n _>oo P(s 7^ s) = 0. 
In the scheme of (T3), the receiver computes 

y = [a H Y-b H CD] mod A 

= [s + ^(HC, b, a)] mod A (6) 

where a G C Mxl and ^(HC.b.a) = (a H H - b H )CU + 
a H Z denotes the effective noise, including the non-integer 
self-interference (due to the fact that cc H H ^ Z[j] lxM in 
general) and the additive Gaussian noise term. The scaling, 
dither and modulo-A operation in d6| is referred to as the 
CoF receiver mapping. Choosing or = b H H _1 , the variance 
of the effective noise is given by cr^, = || (H _1 ) H b|| 2 . This 
scheme is known as exact IFR fTT) . In this way, the non- 
integer penalty of CoF is completely eliminated. More in 
general, the performance can be improved especially at low 
SNR by minimizing the variance of z. rf (H, b, a) with respect 
to a. In this case, we obtain: 

cr 2 (HC,b) = b H C(SNR; 1 I + C H H H HC)- 1 C H b.(7) 

Since a is uniquely determined by HC and b, it will be 
omitted in the following, for the sake of notation simplicity. 




Fig. 2. A deterministic noise-free 2x2x2 finite field interference channel. 

From |T3|, we know that by applying lattice decoding to y 
given in (fSJ the following computation rate is achievable: 



i?(H,b,SNR cfl ) = log + (SNR,„>'(HC,b)) 



(8) 



where log + (a;) = max{log(a;), 0}. Also, the receiver can 
reliably decode M linear combinations S = [BCT] mod A 
with integer coefficient vectors {b^ 1 : £ — 1, ... , M} (i.e., the 
£-th row of B) if 



R < min{i?(H,b £ ,SNR t . fr )} = i?(H, B, SNR olr ) 



(9) 



Using the lattice encoding linearity, the corresponding M 
linear combinations over ¥ q for the messages are obtained 

as 



U = ff -1 ([B] modpZ\j})g- 1 {[C] mod pZ[j])W 

(10) 



= [B],[C],W, 



where we use the notation \B] q = g 1 ([B] mod pL[j]). 

III. Aligned PCoF 

In this section, we propose a novel lattice strategy called 
"Aligned" Precoded CoF (PCoF). This consists of two phases: 
1) The CoF framework transforms a Gaussian network into 
a finite field network. 2) A linear precoding scheme is used 
over finite field to eliminate the end-to-end interferences (see 
Fig. 131. While using the CoF framework, the main perfor- 
mance bottleneck consists of the non-integer penalty, which 
ultimately limits the performance of CoF scheme in the high 
SNR regime (15). To overcome this bottleneck, we employ 
signal alignment in order to create an "aligned" channel matrix 
for which exact integer forcing is possible, as seen in Section 



II-B Namely, we use alignment precoding matrices Vi and 
V 2 at the two sources such that 



[ FfciVa F fc2 V 2 ] 



where H Rfc G 



iMx M 



(ID 



and C Rfc G Z[j] Mx2M . However, the 
precoding over C may produce a power penalty due to the non- 
unitary nature of the alignment matrices, and this can degrade 
the performance at finite SNR. In order to counter this effect, 
we use the concept of IFB (see [12]). The main idea is that 
Vfe can be pre-multiplied (from the oft) by some appropriately 
chosen full-rank integer matrix since its effect can be 
undone by precoding over ¥ q , using [Ak] q . Then, we can 
optimize the integer matrix in order to minimize the power 
penalty. The detailed procedures of Aligned PCoF are given 
in the following sections. 



A. CoF framework based on signal alignment 

In this section we show how to turn any 2-user MIMO IC 
into a noiseless finite field IC using the CoF framework. We 
focus on the first hop of our 2x2x2 network since the same 
scheme is straightforwardly used for the second hop. Consider 
the MIMO IC in Let {yv i e g Fj : < = 1, . . . , M} denote 
the messages of source 1 and {w 2 e eFj : £ = 1, . . . , M— 1} 
denote the messages of source 2. All transmitters make use of 
the same nested lattice codebook C = Ai n Va, where A 
has the second moment a\ = SNR dl . Also, we let Vi = 
[vi i • • • v l<M ] G C x and V 2 = [v 2 ,i • • • v 2 , M -i] € 
C x 1 denote the precoding matrices used at sources 1 
and 2, respectively. They are chosen to satisfy the alignment 
conditions, given by 



F21V1 * 



Fi 2 v 2i fc 
F 22 v 2 1. 



(12) 
(13) 



for k = 1, . . . , M — 1. The feasibility of the above conditions 
was shown in [9] for any M. 

Let Ai G Z[j] MxM and A 2 G J^M-ixM-i denote the 
full rank integer matrices. They will be optimized to minimize 



the power penalty in Section IV Each source k precodes its 
messages over ¥ q as 



= [Afc]- 1 ^, fc = l,2. 



(14) 



Then, the precoded messages are encoded using the nested 
lattice codes. Finally, the channel input sequences are given 
by the rows of: 

X fc =V fc A fe T fc =V,T' fc , (15) 

where t fe t = [/(w' fe £ ) + d fc e ] mod A, and where T k = 
AkT k . Due to the power constraint, each source k must satisfy 

SIMR eff tr(V fc A fe A£v£) < MSNR. (16) 

Since sources use the same nested lattice codes, we can 
choose: 

SNR cfr = min{SNR(V fc ,A fc ) : k= 1,2} (17) 

where SNR(V fc , A fe ) = MSNR/tr(V fc A fc A£V"). 

The decoding procedure is as follows. We first consider the 
aligned received signals. Relay 1 observes 



— Rl = 



FnXi + Fi 2 X 2 + Z Ri 

tl,2 + £2,1 



F u Vi 

— H R , 



H Rl C Rl 



il,Af 



E2.M-1 



(18) 



(19) 



where (a) follows from the fact that the precoding vec- 



tors satisfy the alignment conditions in ( 13 1 and C Rl 
[ Ai C12A2 ] with 

q1xA/-1 
rM-lxM-1 



(20) 



Similarly, relay 2 observes the aligned signals: 



F21V1 

— H Ro 



Hr 2 C R2 



V 4- 1' 

' ^2,1 



t', 



1,M-1 ' i2,M-l 
+' 

5.1, M 



-R2 



Z R , (21) 



(22) 



messages with a change of sign (multiplication by —1 in the 
finite field) at destination 2 (see Fig. [2). Notice that the system 
matrix is fixed and independent of the channel matrices, since 
it is determined only by the alignment conditions. 

Lemma 1: Choosing precoding matrices Mi and M2 as 



Mi = ( jMxM e (_Q 12 Q 21 )) 
Mo - CTM-lxM-l 



-(r 



(-Q 2 iQi 2 )) _ 



(28) 
(29) 



where C R2 = [ A 1 C 22 A 2 
C22 = 



with 



|M-lxM-l 
nlxiW-1 



the end-to-end system matrix becomes a diagonal matrix: 

jMxM q 



(23) 



Q 



Mi 





M 2 



Q = 







1 



M-lxM-l 



The channel matrices in ( fT9] i and ( |22| ) follows the form 
considered in Section |II-B| Following the CoF framework in 
([9} and ([10), if i? < i?(H Rfc C Rfc , B Rfc , SNR dr ), the relay fc can 
decode the M linear combinations with full-rank coefficients 
matrix B„ ,: 



w; 



= [B R J, [ [AJ, [C fe2 ] g [A 2 
( => [B R J 9 [I^ [C M ],] 



wi 



Proof: See the long version of this paper fl6) . ■ 
Based on the above, we proved the following: 
Theorem 1: For the 2 x 2 x 2 MIMO IC defined in {[} 

and (|2), Aligned PCoF can achieve the symmetric sum rate of 

(2M — l)i? with common message rate 

R = min {i?(H Rfc C Rfe , B Rfc , SNR dI ), R{H k C k , B k , SNR'J} 

fc— 1.2 

for any full rank integer matrices A Rfc , B^, B Rfc , and any 



w 9 



matrices Vfe, V Rfc to satisfy the alignment conditions in (13i 
where 

H Rfc = PfciVi, H fe = G fc iV Rl 
SNR d , = min{SNR(V fe , A fc ) : k = 1, 2} 
SNR rff = min{SNR(V Rfc ,A R J : fc = 1,2}, 



where (a) is due to the precoding over ¥ q in ( 14 1. Let W x = 
[BsJ" 1 !!! and W 2 be the first M - 1 rows of [BrJ- 1 ]^. 

Then, we can define the deterministic noiseless finite field IC: and = I Al C fe2A 2 ] and C k = [A Rl C^A, 



Wi 

W, 



= Q 



w 9 



(24) 



where the so-called system matrix is obtained using Ci and 
C 2 as 



Q 



jMxM 
Q21 



Ql2 

jAf-lxM-1 



(25) 



with 

C 12 ,C 22 in ((20) and Q 
Showing that i? grows as logSNR yields: 

Corollary 1: Aligned PCoF achieves the 2M — 1 DoF for 
the 2 x 2 x 2 MIMO IC when all nodes have M multiple 
antennas. 

Proof: See the long version of this paper fl6). ■ 



where 

Q12 = 







lxM-l 



,Q 



21 



rM-lxM-l 



iM-lxl 



jAf-lxM-1 

B. Linear precoding over deterministic network 



IV. Optimization of symmetric sum rates 

Suppose that precoding matrices , V Rfc are determined. 
We need to optimize all integer matrices in Theorem [TJ to 
maximize the sum rates. First, the power-penalty optimization 
problems take on the form: 



Eq. ( 24 1 defines the first hop of the noiseless finite field IC. 
Next, we focus on the second hop. The relay fc uses a precoded 



version of decoded linear combinations W„ 



MfcWt as its 



_ Rfc - -*-k±JL k 

messages. Operating in a similar way as for the first hop, the 
second hop noiseless finite field IC is given by 



argmm 
subject to 



trfVAA H V H l 



A is full rank 



M 

E 



Va, 



(30) 



W R 

w„ 



Q 



w 

-i-l-R 

w„ 



(26) 



where denotes the ^-th column of A. Also, the minimization 
problem of variance of effective noise consists of finding an 
integer matrix B solution of: 



Concatenating ( 24 1 and ((26), the end-to-end finite field noise- 
less network is described by 



argmin max{b£ (S N R~ 1 1 + H H H) " 1 h t } 
subject to B is full rank 



(31) 



W 

-i-l-R 

w 



Q 



Mi 






M 2 



Q 



w, 



(27) 



Lemma [T] shows that the linear combinations decoded at desti- 
nation 1 are equal to its desired messages and are equal to the 



where H denotes an aligned channel matrix and is the ^-th 
row of B. By Cholesky decomposition, there exists a lower 
triangular matrix L such that 

b^SNR" 1 ! + H H H)" 1 b < > = ||L H b £ || 2 . (32) 



^— ^— Time Sharing 

■ — — Aligned PCoF (w/o opt) 

Aligned PCoF 




For the simulation, we used the following precoding matrices 
to satisfy the above conditions: 



V, = [ F^Fial 



F^ 1 1 F 22 1 ] and v 2) i = 1. 



Fig. 3. Performance comparison of Aligned PCoF and time-sharing with 
respect to ergodic symmetric sum rates. 



We notice that problem ( [30] ) (or pT) ) is equivalent to finding 
a reduced basis for the lattice generated by V (or L H ). In 
particular, the reduced basis takes on the form VU where U is 
a unimodular matrix over Z[j]. Hence, choosing A = U yields 
the minimum power-penalty subject to the full rank condition 
in (30 1. In practice we used the (complex) LLL algorithm, 



with refinement of the LLL reduced basis approximation by 
Phost or Schnorr-Euchner lattice search (see the long version 
for details (16)). 

A. Numerical Results: Rayleigh fading 

We evaluate the performance of Aligned PCoF in terms of 
its average achievable sum rates. We computed the ergodic 
sum rates by Monte Carlo averaging with respect to the 
channel realizations with i.i.d. Rayleigh fading ~ CAf(0, 1). 
For comparison, we considered the performance of time- 
sharing where IFR is used for each M x M MIMO IC. 
We used the IFR since it is known to almost achieve the 
performance of joint maximum likelihood receiver [11] and 
has a similar complexity with Aligned PCoF. In this case, an 
achievable symmetric sum rate is obtained as 

R = min{i?(F fcfc , B 1; 2SNR), R(G kk , B 2 , 2SNR)} (33) 

for any full-rank matrices B x and B 2 . The integer matrices 
are optimized in the same manner of Aligned PCoF. Also, 
in this case, we used the 2MSIMR for power-constraint since 
each transmitter is active on every odd (or even) time slot. 
For Aligned PCoF, we need to find precoding matrices for 



satisfying the alignment condition in (13 1. For M = 2, the 
conditions are given by 



FnVi 2 — Fi 2 v 2 i and F 2 iVi i — F 22 v 2 i. 



Also, the same construction method is used for the second 
hop. Since source 1 (or relay 1) transmits one more stream 
than source 2 (or relay 2), the former always requires higher 
transmission power. In order to efficiently satisfy the average 
power-constraint, the role of sources 1 and 2 (equivalently, 
relays 1 and 2) is alternatively reversed in successive time 
slots. In Fig. [3] we observe that Aligned PCoF can have the 
SNR gain about 5 dB by optimizing the integer matrices for 
IFR and IFB, comparing with simply using identity matrices. 
Also, Aligned PCoF provides a higher sum rate than time- 
sharing if SNR > 15 dB, and its gain over time-sharing 
increases with SNR, showing that in this case the DoF result 
matters also at finite SNR. 
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